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Abstract. Patterns of geographic variation in 13 populations of Apodemus agrarius 
from Kaliningrad (GUS) to Macedonia were investigated by means of skull morphology 
(14 variables) and in a subset of 4 populations by electrophoresis (44 enzymes encoded by 
57 gene loci). Genetic distance analysis of biochemical data failed to indicate clusters of 
populations differentiated at the subspecific level. Morphological differences were mainly 
size-dependent. Linear skull dimensions could be attributed to non-genetic, environmental 
adaptations with the exclusion of molars which seem to be relatively invariable against en- 
vironmental conditions. Selective constraints to modify parts of the dentition seem to 
require stronger changes in the genetic program that may vary between different popula- 
tions to a low degree. Looking at all results, А. а. kahmanni shows convergent size relation- 
ships to A. a. istrianus. A. a. kahmanni is in geographic contact with populations of the 
nominal race, and its larger cranial proportions are possibly a result of clinal size variation. 
By contrast, A. а. istrianus is geographically isolated and appears to establish specific 
genetical characteristics as expressed by a highly significantly reduced heterozygosity and 
morphological features similar to those of A. a. kahmanni. 
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Introduction 


The striped field mouse (Apodemus agrarius) inhabits a wide geographical range 
between central Europe in the west and China and Korea in the east (Musser & 
Carleton 1993). In Middle Europe three subspecies of Apodemus agrarius have been 
discussed: A. а. henrici von Lehmann, 1970 from Germany, regarded by some 
authors (e.g. Bóhme 1978) as a synonym of A. а. agrarius, A. a. istrianus Krystufek, 
1985 from Slovenia, and A. a. kahmanni Malec & Storch, 1963 from Macedonia. 
While A. a. kahmanni is regarded as valid by most authors (Bóhme 1978, Kahmann 
& Einlechner 1992), the status of A. a. istrianus was recently questioned by Kahmann 
& Einlechner (1992). 

A. a. henrici was described from Germany (v. Lehmann 1970). Although we had no 
material from the type locality of A. a. agrarius in Russia, we follow Bóhme (1978) 
in synonymizing henrici with agrarius. A. а. istrianus occurs in Slovenia and NE Italy 
(Krystufek 1985, 1991, for Italy see Sala 1974 and Zulian 1987). According to 
Krystufek (1985, 1991, pers. comm. 1995) its populations are geographically 
separated by a gap from east Slovenian populations which represent A. a. agrarius. 
Kahmann (1961) reported on findings from Ribnica, a place right between the two 
current areas, but he left no voucher specimens and Krystufek (1985) could not con- 
firm this locality after intense collecting. Other authors, however, suggested that all 
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A. agrarius from the area of former Yugoslavia and NE Italy should be referred to 
subspecies kahmanni (Djulic & Vidinic 1964, Ondrias 1966, Soldatovic et al. 1971, 
Kahmann & Einlechner 1992). 

Descriptions of subspecific divergence among populations of the striped field 
mouse in Europe were to a great extent based on external morphological traits, 
mainly differences in size. In this paper, we compare patterns of morphological 
differentiation among populations assignable to the 3 subspecies currently recognized 
to their patterns of biochemical differentiation, in order to account for genetic 
relationships that define evolutionary units such as subspecies. Inasmuch, we follow 
the concept of Smith & Patton (1988) to consider those entities to have both 
character (morphological and genetical) and geographic continuity as appropriate 
infraspecific units to be recognized in a formal taxonomy. While from the Oriental 
range of the species only little karyotypic (Bulatova et al. 1991) and biochemical data 
are available (Wang 1985, Zhao & Lu 1986, Liu et al. 1991), the scarce data on Euro- 
pean populations are widely scattered in the literature (Britton-Davidian et al. 1991; 
Filipucci 1992; Gemmeke 1980; Gill et al. 1987; Hartl et al. 1992; Niethammer 
unpubl.). But, dealing with small sample sizes, they seem not to be sufficient to fully 
characterize infraspecific genetic variability of A. agrarius. The purpose of this study 
was to assess the taxonomic status of European populations of A. agrarius at the 
border of its range in western Europe. The present multivariate examination of skull 
proportions in combination with a rigorous analysis of protein variation should give 
answers whether certain population groups warrant recognition as subspecies or not. 


Materials and methods 


Morphometry 

Measurements: In the craniometric part of the study we examined a total of 158 skulls stemm- 
ing from 13 populations between Kaliningrad (GUS) in the north and Lake Dojran 
(Macedonia) in the south (Fig. 1). Only young adult and adult specimens of both sexes (tooth- 
wear classes 3—5 according to Adamczewska-Andrzejewska 1973) were measured in order to 
reduce variance bias in size and shape introduced into the samples by ontogenetically caused 
variation. The sexes were not separated (populations sampled and abbreviation codes are given 
in the legend to Fig. 1). 

Skulls are stored in the following collections: Zoologisches Forschungsinstitut und Museum 
Alexander Koenig, Bonn (ZFMK); Senckenberg Museum, Frankfurt/M. (SMF); Slovene 
Museum of Natural History, Ljubljana (PMS); Staatliches Museum für Naturkunde, Górlitz 
(MNG); private collection H.-J. Pelz, Münster (CP); private collection H. Meinig, Wer- 
ther/Westf. (CHM). 

14 measurements were taken, measurements 1 to 9 (Fig. 2) with a digital calliper (Mitutoyo 
digimatic) to the nearest 0.01 mm, measurements 10 to 14 with a binocular (Zeiss GSZ) with 
an enlargement of 50. All measurements were taken by one of us (H.M.). Abbreviations used 
are: Cbl — condylobasal length (1), zBr — zygomatic breadth (2), IoC — interorbital constric- 
tion (3), RoM — rostral breadth (4), NL — nasalia length (5), MBr — mastoid breadth (6), 
APF — length of anterior palatine foramen (7), MxT — maxillary tooth-row length (8), D 
— diastema (9), MIL — length of first upper molar (10), MIBr — breadth of first upper molar 
(11), M2Br — breadth of second upper molar (12), M3Br — breadth of third upper molar (13), 
ID — incisive diameter (14). 


Statistical analyses | | 
Population genetic measures: Allelic frequencies were computed for each population derived 
from individual electrophoretic genotypes by gene-counting as implemented in the BIOSYS-1 





Fig. 1: Geographical origin of the populations examined: 1 — Zehlau, Kaliningrad area 
(KAL), GUS (7); 2 — Prenzlau, Brandenburg (PRE), Germany (12); 3 — Berlin (ber), 
Germany (13); 4 — Harz, Lower Saxony (har), Germany (15); 5 — Görlitz, Saxony (goer), 
Germany (18), 6 — Osthessen (ohe), Germany (16); 7 — Tiszacsege, Hortobagy (tis), Hungary 
(7); 8 — Radenci, Mura rijeka (rad), Slovenia (21); 9 — Brezice (BRZ), Slovenia (6); 10 — 
Ajdovscina (AJD), Slovenia (11); 11 — Rovinj (rov), Croatia (15); 12 — Banja Bansko (bba), 
Macedonia (3); 13 — Lake Dojran (doj), Macedonia (13). Codes for populations studied 
morphologically and biochemically are given in capitals, codes for populations studied only 
morphologically are given in lower letters; the numbers of skulls measured are given in 
parentheses. 
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program of Swofford & Selander (1981); allele frequency estimates for an isofemale Fl line 
sample from Kaliningrad area (KAL) were corrected for introduced bias not exclusively 
screening for polymorphism in samples from the wild (Long 1993). The amount of genetic 
divergence between populations was computed by Nei's unbiased standard genetic distance D 
(Nei 1978). A phenogram of the genetic relationships among populations was obtained 
performing the unweighted pair group arithmetic average cluster analysis (UPGMA, Sneath 
& Sokal 1973). Standard errors on each bifurcating node were calculated as the standard 
deviation of all pairwise distances between all OTUs joining the nodes within the cluster con- 
secutively (Nei et al. 1985). 

Cranial morphometric analyses: Morphological relationships among geographic samples 
were assessed by four substantial techniques utilizing several statistical routines of the SYSTAT 
version 5.03 for DOS (Wilkinson 1990), the BMDP-PC90 package (Dixon 1990) and the 
NTSYS-pc ver. 1.60 (Rohlf 1990) for IBM-compatible computers. 

Techniques for verification of natural groupings (in this case subspecies) should have the 
property not to be biased by information of group membership, that is an a priori assignment 
of specimens to these groups (Humphries 1984). As an exploratory technique for discovering 
structure in data the Principal Component Analysis (PCA) is widely used in systematic 
studies. Here we employ Multiple Group Principal Component Analysis (MGPCA; Thorpe 
1983, 1988). It provides a multivariate means to assess the within-group components of 
character variation when using intercorrelated linear measurements. By pooling the within- 
group variance-covariance matrices derived from log-transformed cranial variables it contri- 
butes better to among-group discrimination than ordinary PCA. The logarithmic trans- 
formation makes the covariance matrix independent of scaling of measurements but standar- 
dizes variances and preserves allometries (Jolicoeur 1963). Extracted principal components are 
interpreted as patterns of covariation in size and shape, but actually do not confuse the within- 
and between group differences when several groups are used (Thorpe 1976). The first MGPCA 
axis derived from the pooled within-group variance-covariance matrix can be interpreted as 
a general within-group allometric “size” vector if most of the original variables contribute 
with positive signs and equal magnitude to its eigenvector coefficients (Patton & Smith 1990). 

The first step of the procedure was the computation of character residuals from the log- 
transformed variables for each population sample derived from an analysis of variance using 
the MGLH routine of SYSTAT. An ordinary PCA on the covariance matrix of these residuals 
produced eigenvectors to be cross-validated by multiplying the score coefficients with the log- 
transformed variables (using SYSTAT’s weighting variable option). Alternatively, computation 
could be done using BMDP-PC90 tools. First the variance-covariance matrix was computed 
for each of the 13 groups (= populations), and these were pooled to produce a single within- 
group variance-covariance matrix using BMDPAM-module. Then from this matrix the princi- 
pal components were extracted by means of the BMDP4M.-routine. 


The resulting component scores were used in bivariate plots in an attempt to separate the 
groups (= populations or subspecies) either “size” included or excluded (omitting MGPC-1 
= “size-out” analysis). 


Following these latter consideration of a “size-out” analysis (Thorpe et al. 1982), the “size- 
dependent” principal components (MGPC-1 and also MGPC-2) were excluded from sub- 
sequent analyses and the component scores of the MGPCA2-14 res. MGPCA3-14 variates are 
regarded as size-independent ‘characters’, which were subjected as new variables to a discrimi- 
nant analysis to assess grossly size-free variation between populations. Individual scores on 
the first two canonical axes plotted against each other show size-independent shape variation 
among the populations. 


In a slightly different approach used as an independent means to subsume for effects of 
overall size on variation found among populations, cranial variables were first size-adjusted, 
using Burnaby's (1966) canonical variate analysis framework. Data were projected onto the 
hyperplane orthogonal spaced to the *size*loaded vector of the first principal component 
employing the ORTH option of the PROJ module of NTSYS. Individual scores on the 
adjusted principal components plotted against each other show size-independent discrimina- 
tion of the populations. 
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Fig. 2: Skull of Apodemus agrarius with the cranial measurements 1 to 9 indicated (measure- 
ments 10 to 14 not shown). For abbreviations see text. 


Linear Discriminant Function Analysis using the pooled variance-covariance matrix was 
performed to compute the distances between different samples maximizing the between-group 
versus the within-group variance. It requires a beforehand allocation of individual specimens 
to one of the a priori determined groups (Neff & Smith 1979). We graphically demonstrate 
the differences between the groups (= populations) by a Neighbour-Joining tree (cf. Nei 1987) 
clustering the Mahalanobis distances of individual canonical variable scores from group cen- 
troids. Finding classification functions was computationally realized with the ‘Stepwise Discri- 
mination Analysis BMDP-subroutine 7M’. Clustering was done with NTSYS. 

Size and shape covary, and unless isometry pertains, such covariation implies a changing 
relationship between size and shape (Gould 1966). To study this finally, multivariate static 
allometric coefficients for the 14 cranial variables were calculated to look at the influence of 
covariation of shape and form dimensions related to size differentiation (Leamy & Bradley 
1982). 

In a first step we performed principal component analyses separately for each population 
sample (Smith & Patton 1988). Because the first principal component (PCI) of our data satis- 
fies interpretation as a general size factor, the position (= score) of an individual on РС! is 
a measure of its overall body size, while the “raw” loadings (= elements of the eigenvector) 
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of variables on this component describe the relative contribution of each variable to change 
in general size, thus are proportional to allometric coefficients of the characters with respect 
to size (Bookstein et al. 1985). The first principal component of the variance-covariance 
matrix from log-transformed data should therefore represent some kind of an isometric size 
vector that can be rescaled to the length of one (Somers 1986) if covariation between the 
variables approach equality. Where allometry exists, it thus provides a standard measure 
against which growth trajectories of individual cranial characters can be compared (Smith & 
Patton 1988). 

To “normalize” the first principal component to unity we divided its raw loadings by a value 
У, k], where К = raw loadings and p = number of cranial variables) such that their 
squared elements sum up to unity. Then the normalized loadings were divided by l/yp to 
rescale the loadings to be expected if all dimensions (p = 14) have grown at the same rate 
(Shea 1985). Resulting positive allometric variables with multivariate adjusted coefficients > 
1 are those that are relatively larger in large individuals than in smaller ones; negative allo- 
metric variables (coefficients <1) are those with the opposite relationship (Strauss & Book- 
stein 1982). 

Allometric coefficients were used as new variables in a discriminant analysis (employing the 
MGLH routine of SYSTAT) that treats population samples separately. Canonical variable 
plots (Fig. 5) give insight into grouping patterns. 


Electrophoresis 

A total number of 53 animals were caught with snap traps at four localities (no. 1, 2, 9, 10 
in Fig. 1). Tissue samples (muscle, liver, heart) were taken in the field and stored in liquid 
nitrogen until being returned to the ZFMK biochemical laboratory, where they were cut into 
small pieces and maintained in an ultracold freezer (— 85 °C) for long term storage (tissue col- | 
lection). 

Prior to electrophoretic analysis a fivefold volume of 0.1 M Tris/HCl homogenate buffer 
(pH 7.0) containing 0.002M EDTA and 0.05M NADP was added to the weight of portioned 
tissue, either pure organ specific probes or mixes from both liver and muscle, which were then 
homogenized with a motor-driven homogenizer (Polytron dispenser with 12mm shaft, Kine- 
matica, Switzerland) keeping samples cool in an ice-bath. Homogenates were shaken with 0.1 
— 0.2ml Toluene and immediately centrifuged for 10 minutes at 13.000g (Biofuge 13, Her- 
aeus-Sepatech, Germany). The clear supernatant (25/ per sample) was transferred onto Micro 
Test Tissue Culture Plates (COSTAR, Cambridge; Greiner, Germany) and refrozen in a —20 ? 
freezer until electrophoretically processed. 

We employed the procedures of vertical starch gel electrophoresis first described by Smithies 
(1955) and recently reviewed in Geiger (1990), who also gave details due to technical novelties 
and apparative equipments. Starch gels are made in concentration of 12 % and 12.5 % (w/v) 
starch in gel buffer using BIOMOL starch (Hamburg, Germany; Tàb. 1). Handling and prepa- 
ration of gels follows the outlines made by Murphy et al. (1990). Sample application in the 
vertical apparatus is done by means of an Eppendorf comforpette pipetting amounts of 
5—10ul per individual into a preformed slot (20 in total) in the gel, which is then sealed by 
molten vaseline. Gels were electrophoresed overnight (16 h) at 3—4 V/cm in a 4 °C temperated 
freezer, the gels additionly connected to an cooling system with cooling plates. Each gel was 
then sliced into 1.2 mm thick slabs for histochemical overlay-staining adopting the visualiza- 
tion techniques as described by Ayala et al. (1972), Catzeflis et al. (1982), Filipucci et al. (1987), 
Harris & Hopkinson (1978), Hartl & Hóger (1986), Selander et al. (1971) and Shaw & Jain 
(1970). 

44 enzymes and general proteins encoded by 57 presumptive structural gene loci were exami- 
ned for all populations. Electrophoretic running conditions, separation buffer systems used, 
enzymes assayed and their tissue sources are listed in Tab. 1; although no progeny testing was 
routinely done (with the exception of the Kaliningrad area sample KAL) to confirm the mode 
of inheritance of allozyme variants, resulting zymograms generally conformed with simple 
patterns of codominant Mendelian inheritance, so that genetic interpretation of banding pat- 
terns could easily be done based on principles published by Csaikl (1985), Harris & Hopkinson 
(1978), Hartl et al. (1988), Richardson et al. (1986) and Selander et al. (1971). Designation of 
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cumulative coefficient of variation (cv) 





KAL PRE ber har асе ohe tis rad BRZ AJD rov bba doj 


populations 


Fig. 3: Banding diagram of coefficients of variations within populations. 


encoding loci and allelic variation of the allozymes are as follows: Genes are symbolized by 
italicizing the enzyme and protein abbreviation of Table 1; numerical suffixes distinguish 
among multiple zones of cathodal or anodal or both activities on certain zymograms in order 
of decreasing mobility from the most anodal one considering anodal migration first; electro- 
morphs (interpreted as alleles) were given letters in alphabetical order, arbitrarily starting with 
the one that migrated the least to the anode (anodal migration) or the least to the cathode 
(in case of cathodal migrating) under standard electrophoretic conditions as described here 
(Tab. 6). 


General statistical tests 

Modified Mantel’s (1967) randomization test in a multiple regression and correlation extension 
was used to test for matrix associations between genetic, geographical and morphological 
distances among the four populations KAL, AJD, BRZ and PRE, where the distances in one 
matrix are regressed on the distances in the other matrices (Manly 1991). Significance of corre- 
lations between geographic and morphometric distance for all 13 populations in the morpho- 
metric study were tested with ordinary Mantel analysis (1967). 


Results and discussion 
Craniometric analyses 


Variation of single variables 

Coefficients of variation evidence very low intra-populational differences. The 
banding diagram (Fig. 3) shows values as low as 0.018 for Cbl in sample Lake Dojran 
(doj) and a higher value of 0.052 in Osthessen (ohe). As a representative of tooth 
variables M3Br ranges from 0.052 in Kaliningrad (KAL) to 0.099 in Osthessen. The 
diagram shows no disruptive geographical trend due to a characterization of certain 
populations. 


Variation in size 
Condylobasal length (Cbl) and zygomatic breath (zBr) can be considered the most 
useful single indicators of overall cranial size among the variables examined. They 
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are highly correlated with the other skull measurements (less with dentition varia- 
bles; Tab. 2) and have low within-population coefficients of variation (Fig. 3). For 
example, condylobasal length means range from 21.46 mm in population Harz, Ger- 
many (har) to 24.75 mm in population Lake Dojran, Macedonia (doj), representing 
a 13.3 % difference among localities. 

Although our study is faced with a relatively low degree of variability (Tab. 2), 
multiple group principal component analysis was effective enough to discriminate 
between minor morphometrically mensurable differences in cranial size and shape. 

In order to analyze size variation among populations in a multivariate treatment, 
the first two multiple group principal components from the pooled within-group 
character relationships can be considered as general size factors, since all vector coef- 
ficients are positive (tooth variables excluded) and show correlations with the origi- 
nal log-transformed character values (Strauss 1985). The correlation between Cbl, 
for example and MGPC-1 is 0.936. Communalities of the variables that are the pro- 
portions of variance accounted for by the two main factors are given in Tab. 2. Linear 
skull measurements and tooth variables show almost complete loadings on both 
components. 

To investigate the relationships in the craniometric variables on their own, Table 
3 gives the loadings for the three vectors, together with the percentage variation they 
express (cf. Thorpe & Leamy 1982). The first multiple-group principal component 
accounts for 36.48 % of the within-group variation across the entire sampled range 
of A. agrarius in Europe, the first three components account for 69.82 % of total 
variance. МОРС- is the largest (36 o) and is equally loaded in magnitude with con- 
tributions of the cranial variables ID, D, APF, Ком, NL, Cbl and zBr, but inverse 


Table 2: Pearson product-moment correlation coefficients between the mean log;o-transfor- 
med cranial variables for the 13 population samples of A. agrarius, their scores on the first 
three Principal Component axes extracted by a multiple-group PCA and communality of 
variables on the first two components (see text for details). 


character PC-2 communality 


log Cbl 0.176 
log zBR 0-232 
log Ioc 0.457 
log RoM 0.163 
log NL 0.055 
log MBr 0.296 
log APF 0.128 
log MxT 0.505 
log D 0.012 
log MIL 0.846 
log MIBr 0.524 
log M2Br 0.393 
log M3Br –0.003 


eigenvalue 0.001 





% explained variance 23.26 
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correlation by dentition features, emphasizing independence of overall size. The 
second component is dominated by tooth variables (M3Br, M2Br) with less emphasis 
of size contributions from the remaining skull variables. The third component is 
displaying portions of shape variation, expressed through varying and inverse corre- 
lations between variables rendering very low explained variance. Although all varia- 
bles in total contribute most to linear size relationships, there are two principal sets 
of variables, dentition and linear skull measurements, that vary non-concordant 
according to size and shape dimensions. Dentition is mainly independent of indivi- 
dual size in rodents, as expected from the developmental stability of ontogenetic 
growth of the molar dentition, and shows no aptitude to be influenced by environ- 
mental factors. 


Geographic variation of size 

The scores from the first two principal components for the individuals of each popu- 
lation can be used as a multivariate measure of cranial size (Tab. 3). Striped field 
mice have largest skulls in populations Lake Dojran (doj) and Rovinj (rov), and 
smallest in population Harz (har). Mean scores per population increase along the 
MGPC-1 axis from 1.85 in the smallest agrarius populations (har) to ca. 2.15 in the 
kahmanni population (doj) and the istrianus population (rov) (Tab. 6). We found a 
linear population overlap along a mainly size varying array in direction of the 
MGPC.- axis, although the sequence is also mostly influenced by loadings from the 
tooth variables that direct separation along the second MGPC axis (Fig. 4). Means 
of M2Br are smallest in population Osthessen (ohe: 1.05 mm) and largest in Rovinj 
(rov: 1.19 mm), for M3Br the smallest average was again found in Osthessen with 
0.69 mm, and the largest in population Banja Bansko (bba) with 0.83 mm. There is 
a slightly small gap between the agrarius pool and both kahmanni and istrianus sam- 
ples. The same is true for the individual scores grouped according to subspecies 
(Fig. 4). 

In the bivariate plot of the first canonical variable against CV-2 in the “size-out” 
analysis, all differences between the populations are blurred as compared to the size- 
related discrimination described above (Fig. 5). Exactly the same results are gained 
after adjusting the data with Burnaby's discriminant approach, i. e. when most varia- 
tion through size differences between populations were removed from the first princi- 
pal component. Remaining variability left no more clear-cut structure to discrimi- 
nate among groups (not shown here). Eliminating size from the data by means of 
both methods yielded clouds of component scores leaving the populations indistin- 
guishable from one another in character space. 


Discriminant analysis 

A discriminant analysis was performed with populations grouped according to the 
currently recognized subspecies. Canonical variable functions found are useful to 
clearly separate between agrarius, istrianus and kahmanni. The agrarius group is very 
well separated from the istrianus and kahmanni samples, showing small overlap (Fig. 
6). As many as 153 of the 158 animals were correctly assigned to the reference subspe- 
cies (96.2 %). Problems only occurred with the identification of istrianus and kah- 
manni. 
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Fig. 4: Bivariate plot of the 95 % confidence ellipses for the sample means of the first two 
principal component factor scores of 13 populations of A. agrarius derived from a multiple- 
group PCA. The proportion of total variation explained by each component is indicated. The 
inset illustrates character vectors, based on their respective correlations with these axes. 


When the samples were grouped only after their geographic origin, the separation 
of the three subspecies turned out less clearly (Fig. 7); 12 (19 %) animals were ill- 
classified, but 48 % could be classified correctly out of 13 populations, showing a 
high degree of variability. The same result is shown by the Neighbour-Joining tree 
based on the Mahalanobis distances of individual canonical variable scores from 
group centroids. Only a separation of the istrianus/kahmanni group on the one hand 
and agrarius on the other hand can be ensured. The branching pattern within the 
agrarius group displays no significant evidence. 

In general, discriminant analysis is a very useful tool to find or contrast differences 
between groups. However, a pre-allocation to a certain group should always be based 
on hard evidence. It is not advisable to introduce information into the calculations 
that should be confirmed by the following analysis. One should avoid an assignment 
of individuals to subspecies if the purpose of the analysis is to look for subspecific 
differentiation. Results of the discriminant analysis therefore should not be taken as 
affirmative because they are biased by a priori information. In this case multiple 
group principal component analysis revealed a convergent size shaping of W Slove- 
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CV-2 





CV-] 


Fig. 5: Bivariate plot of 95 % confidence ellipses for canonical variate scores on the first 2 
axes derived from a canonical analysis of the cranial “size-out” variables MGPC-3 through 
—14 for the 13 populations of Apodemus agrarius. 


nian and Macedonian populations, not recognized by discriminant analysis. The 
phenomenon of phenotypically similar populations in disparate geographic areas is 
well known (Mayr 1975). 


Static allometric coefficients 

As MGPC-1 accounts primarily for variation in size, its loading reflecting average 
size-related changes among the samples, mean static allometric coefficients likewise 
indicate the manner in which different measurements change in relation to overall 
body size. The allometry values for each of the 14 cranial variables are given in Tab. 
3. We used these allometric coefficients to explore components of genetic programs 
that may underly morphological trait expression and may rule individual growth tra- 
jectories within populations or population groups (= subspecies). To address these 
questions, the static allometric coefficients were used as variables in a discriminant 
analysis with a priori allocation of individuals to the three presumed subspecies. As 
a result, each population sample could be distinguished by a unique set of character 
allometries and their allocation to certain subspecies was with absolute a posteriori 
certainty. However, plotting the first two canonical variables against each other, we 
found three groupings (Fig. 8). Midway lie most agrarius populations, flanked by the 
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Fig. 6: 95 % confidence ellipses surrounding canonical variate scores of the first two discrimi- 
nant functions for the three subspecies groups. 
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Fig. 8: Plot of the 13 populations against their values for two canonical discriminant functions 
derived from the populations! mean static allometric coefficients (see text for details). 
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Fig. 9: UPGMA-phenogram based on Nei’s unbiased genetic distances between 4 populations. 
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istrianus populations which rise to higher values, and the kahmanni populations with 
decreasing values on both axes. The latter group also contains the agrarius popula- 
tion Tiszacsege (tis) from Hungary. 

We again found highly informative allometric relationships to characterize indivi- 
dual populations. Reasons for this may be that A. agrarius can quickly adapt its 
growth to changing environmental conditions. Skull dimensions appear to be very 
easily transformed when different food ressources are exploited, for instance in 
urban green areas (Sikorski 1982). 


Electrophoretic analysis 


Genetic distances (Nei's unbiased standard measure D), sample size per locus, 
percentage polymorphic loci, and direct-count and expected mean heterozygosity 
under a Hardy-Weinberg equilibrium are listed in Table 4. Electrophoresis indicates 
a low level of variation both within and among populations. Of 57 loci analyzed, 
twenty-seven were fixed for the same allele in all populations screened for protein 
variation. The remaining variable loci were each polymorphic for two alleles recom- 
bined in different genotype frequencies (Tab 6). A third allele (c) of the enzyme Gpt 
(rf = 38 mm) was detected in all populations but Ajdovscina. On average, we found 
1.2 alleles per population. 

The number of loci expressing variation within populations ranged from 15.8 9/o 
(KAL) to 22.8 % (BRZ), using a 5 % frequency cutoff level. Direct-count heterozy- 
gosity per population ranged from 0.057 + 0.017 in population Ajdovscina to 0.095 
+ 0.027 in population Prenzlau. All populations display strikingly low values of 
genetic distances (D ranges from 0.009 to 0.040). 

Values of genetic variation are within the range reported in previous work on 
A. agrarius (Filipucci 1992; Britton-Davidian et al. 1991). The latter authors, who 
studied kahmanni populations from Greece and Bulgaria, especially considered 
genetic distances as falling within the values generally recorded for subspecific gene- 
tic differentiation, whereas Filipucci (1991) contrasted this opinion by stating that 
“a relatively low value of genetic distance (D = 0.027) was observed among the 
populations of A. agrarius, which are attributed to different subspecies: A. a. istria- 
nus and A. а. agrarius ..:. 

We conclude from our findings that genetic variability values clearly demonstrate 
an amount of genetic differentiation to be found in local populations of a species 
with a high level of gene flow among conspecific populations. Clusters do not indi- 
cate any significant branches among the populations (see error bars in Fig. 9), as 
must be expected when certain subspecies are involved. A very limited degree of dif- 
ferentiation over its range can often be found in species with distribution patterns 
of a typical Euro-Siberian faunal element (de Lattin 1967), as A. agrarius can be 
described. Zhao & Lu (1986), for example, reported on a similar mobility variation 
in serum proteins among Chinese populations of A. agrarius. Direct count heterozy- 
gosity of population Ajdovscina, however, provides evidence of a genetic disconti- 
nuity concerning a highly significantly lowered mean level of gene diversity (one-tai- 
led tio; 15 = 4.89), as compared to the E Slovenian population Brezice, indicating ап 
isolated gene pool. 
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Fig. 10: The samples AJD, rov, BRZ and rad and the local distribution patterns of A. agrarius 
in the northern Balkans (shaded) (after Krystufek 1991 and Petrov 1992). 


Mantel test on distance matrices 


The relationship between populations of A. agrarius in Europe based on morpholo- 
gical, genetical and geographical distances was studied in a subset of four popu- 
lations. The Mantel test was performed to test for statistical association between 
three distance matrices. Morphological distance is represented by Mahalanobis D? 
between the 4 populations previously derived from the discriminant analysis of the 
log-transformed cranial variables. Nei's unbiased genetic standard distance reflects 
genetic differences, and geographic distances were measured as straight air line 
distances between localities. 

We did not find any high Pearson correlation coefficients to indicate significant 
intercorrelated associations between matrices. The reason may be the restricted data 
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Fig. 11: Scatter diagram of CbL against IoC for W (AJD, rov) and E Slovenian (BRZ, rad) 
populations of A. agrarius. 


set of only four populations, including the Slovenian populations BRZ and AJD that 
show a relatively higher proportion of morphological difference in relation to their 
actual geographic distances (about 120 km), as compared to the other populations 
studied (Fig. 10). When CbL is plotted against IoC, W and E Slovenian populations 
are clearly different (Fig. 11); their morphological divergence is as great as that 
between German and Macedonian populations (Fig. 7). 

Ordinary Mantel test on significance probability of matrix association between 
morphological and genetic distances between all populations gave over 95 9/o correla- 
tion (p = 0.0265), testifying that the the morphological differentiation increases with 
distance between populations. 


Conclusions: body size variation, genetics and systematics 


In phenetic analyses of geographic variation, one looks for geographic character 
patterns which components reflect simple plastic responses to local environmental 
conditions to clearly distinguish them from fundamental adaptive genetic changes 
which fit the requirements of the subspecies concept. In this context a first step is 
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Table 5: Means, standard deviations, coefficients of variation and range of the craniometric 
measurements within 13 populations of A. agrarius. 
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Table 5 (continued) 
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Table 5 (continued) 
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Table 5 (continued) 
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KAL: Zehlau, Kaliningrad area, GUS (п = 7); PRE: Prenzlau, Brandenburg, Germany (п 
= 12); ber: Berlin, Germany (n = 13); har: Harz, Lower Saxony, Germany (n = 15); goer: 
Görlitz, Saxony, Germany (п = 18); оће: Osthessen, Germany (n = 16); tis: Tiszacsege, 
Hortobagy, Hungary (n = 7); rad: Radenci, Mura rijeka, Slovenia (n = 21); BRZ: Brezice, 
Slovenia (n = 6); AJD: Ajdovscina, Slovenia (n = 11); rov: Rovinj, Croatia (n = 15); bba: 
Banja Bansko, Macedonia (n = 3); doj: Lake Dojran, Macedonia (n = 13). 
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Table 6: List of loci, electrophoretic mobility (rf-values [mm]) and interpreted genotypes 
detected in the 4 populations of A. agrarius. 


Prenzlau Brezice Kaliningrad Ajdovscina 
genotypes mobility [rf] (n = 16) (п = 15) (п = 8) (п = 14) 
observed numbers 
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Table 6 (continued). 


Prenzlau Brezice Kaliningrad Ајаоуѕсіпа 
genotypes mobility [rf] (n = 16) (n = 15) (n = 8) (n = 14) 
observed numbers 
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to part morphometrical variation into size and shape components. Size is more likely 
to be affected by fluctuations of the external environment, whereas differences in 
body proportions generally provide more reliable indications of internal, genetically 
controlled, shape building processes (Boone et al. 1993). 

To establish subspecies as infraspecific evolutionary units one ought to find gene- 
tic divergence. From the genetical point of view colonization events in historical 
times can be invoked to explain the low degree of genetic differentiation. A. agrarius 
is supposed to have expanded its range into western Europe from eastern central sett- 
lements since about 7000 A. D. (Bóhme 1978). Time to diverge in the newly occupa- 
ted areas has therefore been too short to generate genetic variation that can be detec- 
ted with genetic distance measures. Gene flow across the populations has probably 
never been interrupted so long. In this study we found significantly reduced hetero- 
zygosity within the W Slovenian population (AJD) as an indication of recent iso- 
lation and the potential for evolutionary independence. 

Looking at heritable portions within morphological traits some features of the 
dentition (M2Br, M3Br) seem to be relatively stable against modificatory adaptabi- 
lity of skull dimensions. Molars have equal size in kahmanni and istrianus popula- 
tions, and these population groups which have no geographical contact are possibly 
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expressing similar genetical characteristics that must have been developed indepen- 
dently in both gene pools. 

The most obvious fact from our study is a clinal size variation that increases from 
north to south with W Slovenian and Macedonian populations showing about the 
same size. These populations have been named as subspecies istrianus (Krystufek 
1985) and kahmanni (Malec & Storch 1963). We actually cannot follow the argumen- 
tation of Kahmann & Einlechner (1992), based solely on size criteria, to synonymize 
the subspecies istrianus with kahmanni, because there exists a distributional gap in 
Istria (Fig. 10) documented by Krystufek (1985, 1991), but a continuous distribution 
from E Slovenia to Macedonia (Petrov 1992). From our point of view, kahmanni sim- 
ply could be the final chain-element of a clinal size variation that suffers from non- 
genetical environmental impact (see Fig. 8: similar allometric growth of Hungarian 
and Macedonian populations). Considering istrianus, there is enough geographical, 
morphological and genetic divergence to warrant subspecific nomination of the NE 
Italian and W Slovenian populations as A. a. istrianus. 

To test the hypothesis that kKahmanni is part of the agrarius gene pool and that 
istrianus is an isolated population group with subspecific status, it may be useful to 
investigate craniometrically and electrophoretically samples originating from the 
area between E Slovenia and Macedonia. 
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Zusammenfassung 


An Stichproben von 13 europáischen Populationen der Brandmaus aus Deutschland, GUS, 
Slovenien, Kroatien, Mazedonien und Ungarn wurden 14 Schädelmaße für eine multivariate 
morphometrische Analyse (Mehrfach-Gruppen-Hauptkomponentenanalyse, lineare Diskrimi- 
nanzanalyse, statische Allometriekoeffizienten) genutzt, um morphologische Differenzierun- 
gen aufzuzeigen, auf deren Grundlage die bisherige infraspezifische Gliederung der Art diku- 
tiert wird. Außerdem wurden 4 Populationen enzymelektrophoretisch untersucht (Berechnung 
genetischer Abstánde auf der Basis von 57 Enzymloci), um ein genetisches Korrelat zur mor- 
phologischen Variabilität zu bekommen, mit dem das Ausmaß des Genflusses zwischen den 
Populationen abzuschátzen ist. Anhand dieser enzymphánotypisch bzw. über Proportionali- 
tätsänderungen der kraniometrischen Variablen aufzeigbaren Unterschiede wird die Nützlich- 
keit subspezifischer Abgrenzungen unter dem Aspekt evolutiver Eigenentwicklungen infolge 
geographischer Isolation bzw. unterschiedlicher Besiedlungsfolgen untersucht. A. agrarius 
weist eine nur geringe genetische Variabilität auf, die kaum Rückschlüsse auf subspezifische 
Differenzierungen zuläßt. Die morphologische Analyse zeigte vor allem größenabhängige 
Differenzierungen mit vermutlich modifikatorisch bedingten Ausprägungen, welche als 
Grundlage für die bisherige subspezifische Gliederung dienten. Alle festgestellten Einzel- 
befunde deuten darauf hin, daß A. agrarius kahmanni aus Mazedonien Kontakt zu Populatio- 
nen der Nominatform hat und wahrscheinlich lediglich das Endglied einer klinalen Größen- 
zunahme in nord-südlicher Richtung darstellt. Dagegen ist A. agrarius istrianus aus Slovenien 
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von dem restlichen Verbreitungsgebiet der Art getrennt, was in einem hochsignifikant herab- 
gesetzten Heterozygotiegrad und einer deutlich morphologischen Differenzierung, konvergent 
zu A. a. kahmanni, zum Ausdruck kommt. 
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